Thermodynamic modelling of phase separation in manganites 
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We present a phenomenological model based on the thermodynamics of the phase separated state 
of manganites, accounting for its static and dynamic properties. Through calorimetric measurements 
on Lao.225Pro.4oCao.375Mn03 the low temperature free energies of the coexisting ferromagnetic and 
charge ordered phases are evaluated. The phase separated state is modeled by free energy densities 
uniformly spread over the sample volume. The calculations contemplate the out of equilibrium 
features of the coexisting phase regime, to allow a comparison between magnetic measurements 
and the predictions of the model. A phase diagram including the static and dynamic properties 
of the system is constructed, showing the existence of blocked and unblocked regimes which are 
characteristics of the phase separated state in manganites. 

PACS numbers: 75.30.Kz, 72.80.Ga, 71.23.Cq 



I. INTRODUCTION 



The issue of phase separation is currently one of the 
main topics of research in strongly correlated electron 
systems, [lj Phase separation (PS) fully develops in vari- 
ous manganese binary oxides, but there are also evidences 
of the key role played by clustered states in high Tc 
superconductors. [2 Yet, after several years of intense ex- 
perimental and theoretical research in this area, the true 
nature of the phase separated state observed in mangan- 
ites is still controversial. Some phenomenological mod- 
els points to the strain between the coexisting phases 
as the main reason for the appearance of phase sepa- 
rated states. 0,0] In addition, there is a lot of theoreti- 
cal evidence in favor of the role of intrinsic disorder in 
the stabilization of the phase separated state. These 
models are based mainly in the double exchange the- 
ory, with a fundamental role played by electron-phonon 
coupling, la, Ifj, [Tj Khomskii and co workers have also 
pointed out the tendency to PS of double exchange 
Hamiltonians when elastic interactions are included. |8j 

The presence of quenched disorder can lead to a rough 
landscape of the free energy densities, triggering the for- 
mation of clustered states, which are induced by phase 
competition. Moreo et al. obtained phase separated 
states in a Monte-Carlo simulation of a random-field Ising 
model when disorder is included in the coupling and ex- 
change interactions. Similar results were obtained by 
Burgy and coworkers, using a uniform distribution of ex- 
change interactions. |lO| As shown for the case of first 
order transitions, quenched impurities can lead, under 
certain circumstances, to a spread of local transition tem- 
peratures (where local means over length scales of the or- 
der of the correlation length) leading to the appearance 



of a clustered state with the consequent rounding of the 
first order transition. Ll| Typical mechanisms to include 
disorder are chemical|l2j and structural. [Til llij 

Despite the intense effort towards a microscopic un- 
derstanding of the manganites. |l5lH^ | many macroscopic 
features of the phase separated state, including its ther- 
modynamic properties and dynamic behavior, still re- 
mains to be studied in greater detail. One of the most 
interesting features of the phase separated state is the en- 
twining between its dynamic and static properties. Some 
of the phase separated manganites display slow relax- 
ation features that hide experimentally the real equilib- 
rium thermodynamic state of the system. |l7| This is why 
the construction of phase diagrams is currently focused 
on the dynamic properties of the phase separated state, 
with regions of the phase diagrams nominated as "frozen" 
or "dynamic" PS,0 "strain glass" or "strain liquid" [lflj 
or, directly, ascribed as spin glass phases. (2(1 l21| 

Among the challenging issues that have not yet been 
addressed is an understanding of the phase separated 
state in terms of its thermodynamic properties. An anal- 
ysis of the behavior of phase separated systems based on 
the probable free energies functional has been schemat- 
ically realized, jlM E3 but without the corresponding 
measurements supporting the proposed scenario. In the 
present study we perform an attempt to construct the 
thermodynamic potentials of the FM and non-FM phases 
of a PS manganite, through calorimetric and magnetic 
measurements. The experiments were carried out in a 
polycrystalline sample of Las/g-j/PryCas/gMnOa (y — 
0.4), a prototypical phase separated system in which the 
effects of the substitution of La by Pr produces an over- 
whelming effect on its physical properties. 0] We took 
advantage of the fact that in the mentioned compound 
homogeneous phases can be obtained at low tempera- 
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tures in long time metastable states, which allow us to 
measure separately the specific heat of each phase, CO 
or FM, in the low temperature region (between 2K and 
60K). With this data it is possible to write an expres- 
sion for the difference between the Gibbs energies of the 
homogeneous phases as a function of temperature and 
magnetic field. The needed constants to link the ther- 
modynamic potentials of both phases were obtained from 
indirect measurements based on the static and dynamic 
behavior of the system. The phase separated state is 
modeled through the hypothesis that the free energy den- 
sities are spread over the sample volume, and that its 
non-equilibrium features are governed by a hierarchical 
cooperative dynamic. Within this framework it is pos- 
sible to construct a phenomenological expression for the 
free energy of the phase separated state based on exper- 
imental data, which is able to describe consistently the 
behavior of the system as a function of temperature and 
applied magnetic field. 



II. EXPERIMENTAL DETAILS 

The measurements were made on a polycrystalline 
sample of Lao.225Pro.4oCao.37sMn03. Details of material 
preparation were previously published. 18] Both magne- 
tization and specific heat results were obtained with a 
Quantum Design PPMS system. Magnetization data was 
measured with an extraction magnetometer, as a function 
of temperature, applied magnetic field, and elapsed time. 
All temperature dependent data was measured with a 
cooling and warming rate of 0.8 K/min. Specific heat 
data was measured with a relaxation method, between 2 
and 60 K. 



III. RESULTS AND DISCUSSION 



Following the H — T phase diagram of the 

compound. |l8l flf| after zero field cooling (ZFC) the 
sample reaches the low temperature state mainly in the 
CO phase, a state that has been described as "frozen 
PS" 0| or "strained glass" This frozen state can 
be released by the application of a moderate magnetic 
field (i?=2.2T),[24j above which the compound trans- 
forms into the FM phase in an abrupt metamagnetic 
transition. [2^ I2H. EtI ] After this step transition the sam- 
ple remains in this homogeneous FM state until a tem- 
perature around 70 K, even after the magnetic field is 
removed. These facts were used to perform the measure- 
ments of the specific heat c a of each phase (a =CO or 
FM) between 2K and 60K, each one considered as ho- 
mogeneous in this temperature range under specific field 
conditions. In Fig. 1 the plot of c/T vs T 2 is shown 
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FIG. 1: (color online). Specific heat of 

Lao.225Pro.4oCao.375Mn03 as a function of temperature 
for H=0 in the CO phase (squares) after ZFC, and in the 
FM phase (triangles) after ZFC to 2K following of a 0-9T-0 
field sweep. The inset shows an enlargement of the low 
temperature data. 



for measurements performed while warming after ZFC 
with different procedures: under zero field (CO phase), 
and under different fields H after a field sweep 0-9T-H, 
for H=0, 1 and 2 T (FM phase). The data of the CO 
phase and that of the FM are clearly distinguished. Also, 
the results obtained for the FM phase for the fields em- 
ployed are practically identical, which is a signature that 
the FM phase obtained after the application of 9T re- 
mains homogeneous until the highest temperature inves- 
tigated. Besides, the fact that the specific heats of the 
FM phase are almost independent of H is a signal that, 
in the range of temperature investigated, there is no sig- 
nificant field dependent contribution to the entropy of 
the FM phase, indicating that the magnetization is satu- 
rated for all fields. The data obtained was adjusted using 
standard models for CO and FM phases. [2S| The small 
upturn observed at low temperatures corresponds to the 
onset of the ordering of the magnetic momentum of the 
Pr atoms. 29] 

The thermodynamic Gibbs potential g of each phase 
may be written as g%{T,H) = E a (T, H) - TS a (T, H) 
where the superscript a indicates the phase (CO or FM), 
E and S are the enthalpy and the entropy respectively. 
From the specific heat data we could construct both E 
and S in the usual way: 123 



E a (T,H) = E c n 



c a (T)dT - ME 



To 



S a (T,H) = 5 Q 



c a (T) 
T 



dT 



(1) 



(2) 



T 
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where M is the magnetization of the phase (we assume 
M=0 for the CO phase) and To=2K is the lowest tem- 
perature reached in the measurements. As the specific 
heat of the FM phase was found almost independent of 
E , we consider the Zeeman term in Eq. (1) as the only 
dependence of the Gibbs potentials of the FM phase with 
E, so that go M (T, E) = g<f M (T,0) - ME. No depen- 
dence with E is considered for the CO phase, since there 
is no way to perform the measurements on the CO phase 
under an applied E due to its instability against the ap- 
plication of E in the temperature upturn. The terms Eq 
and Sq are respectively the values of the enthalpy and 
the entropy at the initial temperature To. Since we are 
interested in the energy difference between the phases 
involved, we have taken Eq M — and Sq M = as refer- 
ence values, remaining Eq° and SQ°as the constant to 
be determined. 

In order to determine Eq° we followed the previously 
published experimental data in relation to the abrupt 
magnetic transition from the CO to the FM phase, which 
happens at low temperatures (below 6K), under a mag- 
netic field of around 2.2 T [Ref. 0]. This transition is 
accompanied by a sudden increase of the sample temper- 
ature, which reaches a value around 30K after the tran- 
sition. Due to the velocity of the process, it is plausible 
to consider that the enthalpy remains constant at the 
transition point, so that the following relation is fulfilled: 

30if 

E£°= J c FM (T)dT — M sat E (3) 

where M sat is the magnetization of saturation of the FM 
phase (M sat =3.67nB /Mn=20.5 J/molT) and #=2.2T. 
This calculation yields Eq° = 28.3J/moZ, indicating 
that the homogeneous FM phase has lower free energy 
at low temperatures even for E—0, as suggested in Ref. 

m 

In order to fully construct the thermodynamic po- 
tentials we have to determine the remaining constant 
Sq° which, at this point, is what controls the transi- 
tion temperature between the homogeneous phases, pro- 
viding it is a positive quantity as expected for the dif- 
ference of entropy between the FM and CO phases due 
to the excess of configurational entropy of the latter. [22^ 
In Fig. 2 the obtained thermodynamic potentials of 
the homogeneous phases are displayed, assuming a value 
Sq° = 0.65J/molK, this value was obtained by adjust- 
ing the M{E) curve of Fig. 4b, as explained below. The 
plot indicates that the homogeneous FM state has a lower 
energy than the CO state for all fields at low tempera- 
tures, while the CO phase is the stable state at high tem- 
peratures, with field-dependent transition temperatures 
ranging from 30K for E=0 to 60K for E=2T. 




■ H(T) ■ 

20 40 60 

T(K) 

FIG. 2: (color online). Free energies as a function of T of the 
homogeneous phases calculated through Eqs. (1) and (2) for 
different field sweeps. Inset: temperature of the sample as 
a function of H when the step transition occurs (taken from 
Ref. 13); the base temperature is 2.5 K. 



The above presented results give us an insight on the 
behavior of the system under the hypothesis that no 
phase separation occurs, i.e., describes the thermody- 
namics of the homogeneous equilibrium phases. In ad- 
dition, in order to obtain a phase separated state from 
the thermodynamic data an appropriate modelling needs 
to include a priori the existence of the phase separated 
state. However, one needs to be careful when comparing 
the predictions of the model with experimental results. It 
is well established |lJ3, ll9| that the phase separated state 
is characterized by a slow dynamics, which implies that 
equilibrium is hardly reached in laboratory times. The 
equilibrium properties must be linked with the measured 
data and therefore a dynamic treatment is needed. 

In the discussion that follows we perform a qualita- 
tive analysis within a framework where both static and 
dynamic properties are treated on a phenomenological 
basis. It is well known that the physical properties of the 
LPCMO system changes dramatically near y=0.32, 0] 
revealing the extreme sensitivity of the system to small 
variations in the mean atomic radius of the perovskite A- 
site. This can be due to the effects that quenched disor- 
der introduced by chemical substitution has on the local 
properties, [|J or else to the role played by " martensitic- 
type" accommodation strains originated by the volume 
differences between the FM and CO unit cells. 0, H3| 
On one hand, the inclusion of disorder in a random-field 
Ising model leads to a spatially inhomogeneous transition 
temperature, from the paramagnetic disordered phase to 
the " ordered" phase, characterized by the appearance of 
clustered states. [lfJ,|3J| This fact implies a spatial depen- 
dence of the quadratic coefficient in a Landau-type ex- 
pansion. On the other hand, strain induced by the shape- 
constrained transformation between the CO and the FM 
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phases could lead to the phase separated state through 
the frustration of long-range interactions. [I3I.I22]] In the 
latter view, the properties of a specific compound are 
governed by an "effective" Pr concentration, which con- 
trols the capability of the system to accommodate the 
anisotropic strain. [22| These two alternative pictures are 
not mutually exclusive; the clustered states induced by 
disorder are enhanced if correlated disorder is included 
in the model, [^3 a feature that mimics the cooperative 
effects of the Jahn- Teller distortion, in a similar way as 
elastic interactions are able to induce long scale phase 
separation in phenomenological models p| (in this last 
case a renormalized fourth order term is the responsible 
for the introduction of spatial inhomogeneities). Addi- 
tionally, local variations of the atomic composition can 
couple with the anisotropic strain triggering the forma- 
tion of the phase separated state. j2^| 

These facts can be qualitatively described introducing 
non-uniform free energy densities for the CO and FM 
phases. The simplest form is a uniform distribution of 
these densities over the sample volume, with mean val- 
ues equal to the free energy of the homogeneous phases. 
Within the hypothesis that precursor effects of phase sep- 
aration are due to variations of the local composition, we 
follow Imry and Wortis theory [llj describing the effects 
of disorder on a system displaying a first order transi- 
tion in order to estimate the width of the free energies 
distributions. Following their ideas, and considering that 
disorder affects mainly the free energy density of the FM 
phase, (T r n is nearly constant as a function of Pr content 
y [Rcf. [12j]) we can write an expression for the local free 
energy density g FM (Ay) depending on the fluctuations 
of composition Ay. 



FM I A \ FM i / c 

9 {Ay) = g + (S 



FM 



-•CO 



iFM 



)f A„ (4, 



where Ay is taken over length scales comparables with 
the correlation length, which is around 1 nm for micro- 
scopic clusters in Prn7Can3Mn03 1331 With the Gaus- 
sian distribution proposed in Ref |llj Ay can be as high 
as 0.1 over nanometer length scales; the development of 
micrometer sized domains would need of the considera- 
tion of elastic interactions. For the sake of simplicity we 
take a uniform distribution for Ay over the sample vol- 
ume. With these assumptions, the free energy density of 
the FM phase as a function of x (the volume coordinate 
is normalized to 1) is written as: 



9* M (x)=9& 



FM 



6(1/2 



(5) 



where 6 is the parameter controlling the width of the 
free energy functional, which can be estimated from Eq. 
(4) as S « (S co - S FM )^y/2 « (70K)SS°, taking 
into account that ^j 2 - « — 350K [Ref.[l3|] and assum- 



ing Ay w |(l/2 — x). In this way, the equilibrium FM 
fraction x eq at the given T and H is obtained as: 



.</,> 



c °(T)-g™(T,H) 



(6) 



This expression could be used for determining the pa- 
rameters S co and 6, using for instance the M(H) data 
at different temperatures. However, as stated before, 
the global behavior of the system at low temperatures is 
characterized by out-of-equilibrium features, so the true 
values for x eq are not easily accessible experimentally. 

In order to circumvent the fact that the thermody- 
namic equilibrium state is not reached experimentally, an 
alternative approach is to consider that the response of 
the system within the phase separated regime as a func- 
tion of temperature can be qualitatively described within 
a model of cooperative hierarchical dynamics, using an 
activated functional form with state-dependent energy 
barriers. The time evolution of the FM fraction x is 
given by: 



dx (| X g q X^ 
dt \Xrn — x\ 



v e 



(7) 



where vq represents a fixed relaxation rate and U is a 
(field dependent) energy barrier scale. This model is 
similar to that employed to describe vortex dynamics 
in high Tq superconductors |34| and is based in dynamic 
scaling for systems with logarithmic relaxations. [35j The 
inbreeding between the dynamic behavior and the equi- 
librium state of the system is given by the functional 



form of the effective energy barriers 



U(H) 



which di- 



\x eq — x\ ' 

verge as the system approaches equilibrium. This fact 
represents the main difference with respect to a pure 
superparamagnetic behavior, and predicts the existence 
of state-depending blocking temperatures Tb(H,xb) at 
which, for any given FM fraction xb , the system be- 
comes blocked, in the sense that the velocity of change of 
the FM fraction is lower than, for instance, the detectable 
velocity ^ cxp , estimated as « 10 _6 s _1 for a conventional 
data measurement that takes 30 sec. 18] . This gives the 
following relation for Tb{H,xb)'- 



T b (H,x b ) < 



U(H) 



| ln(l> /Wexp) 



(8) 



Through Eq. (8) it is possible to obtain an experi- 
mental estimation of the factor U (H) / \yi(vq / v exp ) which 
governs the interplay between the dynamics of the sys- 
tem and the measurement procedure. Figure 3 shows 
M(H) data for selected temperatures, and the experi- 
mentally obtained values of Tb{H 1 xb = 5) as a function 
of H, determined through low temperatures measure- 
ments between 6 and 32 K. The values of H are those 
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FIG. 3: (color online). M(H) data at selected temperatures: 
2.5, 8, 15, 25, and 30 K. Inset: log-log plot of T B versus H/T B , 
where H is that for which M is half of its saturation value. 



at which the system reaches the state with ig = 0.5. 
The main assumption is that each point of the M(H) 
curve corresponds to a blocked state compatible with 
the measurement procedure. As indicated in the inset 
of Figure 3 the relation \T B {H, \) « 25 J ff"' 3 holds, with 
(3 w 2.33. We will show later that the equilibrium state 
at low temperatures is fully FM, so Eq. (8) is a di- 
rect measurement of the field dependence of the factor 
U(H)/ln{v a /v cxp ) = \T B (H,\). Within the dynamic 
model adopted, this factor is independent of the partic- 
ular value of Xb choosen for its determination. With the 
assumption that this relation holds in the whole temper- 
ature range where phase separation occurs, we can write 
a simple equation relating the equilibrium state with the 
experimentally accessible state 



,(T,H)-x B (T,H)\ < 



U(H)/ln{v /v cxp ) 
T 



(9) 



Through this equation it is possible to make the link be- 
tween the equilibrium state of the system, which can be 
obtained through the free energies functional, and the ex- 
perimental data obtained through both M(H) and M(T) 
measurements. 



The upper panel of Fig. 4 sketches the tempera- 
ture evolution of x eq for different magnetic fields, ob- 
tained from Eq. (6) with Sq° = 0.65J/molK and 
S w (70K) SH° = 45.5J/mol. It shows that, with the 
set of parameters employed, the low temperature state 
of the system is fully FM for moderate fields. However, 
the accessible FM fraction after a ZFC procedure is small 
for H < 2 T, due to the weight of the blocking term in 
Eq. (9). This is why, besides the fact that the difference 
between the free energies of the CO and the FM states 
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FIG. 4: (color online), a) Temperature dependence of the 
equilibrium FM fraction, x eq , for the indicated fields, b) FM 
fraction x as a function of H obtained through M(H) data 
(open symbols), and calculated with Eq. (9) (solid line) for 
T — 40K. The FM equilibrium fraction x eq obtained from 
Eq. (6) is also displayed (dashed line), c) Field needed to 
make the system half FM, as a function of T, from M(H) 
measurements (open symbols) and calculated through Eq. (9) 
(solid line). The temperature dependence of the field needed 
to make the equilibrium state of the system fully FM is also 
displayed (dashed line). 



increases as temperature is lowered, the magnetic field 
needed to induce the CO-FM transition also increases, 
a fact that at first sight could be interpreted as a reen- 
trance of the CO state. As the temperature is raised 
above 25K the influence of the blocking term decreases; 
in this temperature region the main factor determining 
the field needed to induce de CO-FM transition lies in 
the field dependence of the equilibrium fraction. 



The middle and lower panels of Figure 4 show the com- 
parison of experimental data with the results obtained 
through Eq. (9). In Fig. 4b, M(H) data at 40 K normal- 
ized by its saturation value and the corresponding calcu- 



lated values are displayed, showing the good agreement 
between them. Figure 4c shows measured and calculated 
values of the field H c at which xb = 0.5, as a function 
of temperature; this is the field needed to make half the 
system FM. Also shown is the temperature dependence 
of the field needed to make the equilibrium state of the 
system fully FM. As can be seen, the calculated curve for 
H c reproduces the experimental behavior, with a mini- 
mum around 30K. This minimum signals the crossover 
from the blocked regime at low temperatures (frozen PS) 
to the coexistence regime at higher temperatures (dy- 
namic PS). In the frozen PS regime the stable state of 
the system is homogeneous FM for all fields needed to 
induce the growth of the FM phase to a value x = 1/2; 
the presence of the CO phase is only explained by the 
slow growing of the FM phase against the unstable CO 
due to the energy barriers. In the dynamic PS regime, 
the influence of the blocking term diminished, the equi- 
librium state is true PS for moderate fields, and the effect 
of H is mainly to unbalance the amount of the coexisting 
phases. 

Figure 5 shows the x — T state diagram obtained from 
Eqs. (6) and (9) for a field H = 1.3T, which displays 
both dynamic and static properties of the model system. 
A line in the phase diagram divides it in two major re- 
gions, depending on the equilibrium FM fraction x eq . A 
point above the line indicates that the system has an 
excess of FM phase; below the line the state of the sys- 
tem is characterized by the presence of metastable CO 
regions. Each of these regions is in turn divided in two 
other, labeled dynamic CO or dynamic FM, indicating 
that if the system is in a state within this region is able to 
evolve toward equilibrium within the measurement time. 
The regions labeled as frozen indicate that the system 
is blocked, and no evolution is expected within the time 
window of the experiment. Data points x(T, 1.3T) ob- 
tained from M(H) and M(T) measurements are also 
showed. The M(H) data is obtained in a field sweep 
after ZFC to the target temperature, and the M(T) in a 
field warming run after ZFC to 2K. The data extracted 
from the M(H) measurements gives information on the 
x values for which the system becomes blocked at each 
temperature, for the specific field and for the character- 
istic measuring time, indicating the frontier between the 
dynamic and frozen FM regions. The data obtained from 
M(T) measurements coincides with that of M(H) in the 
low temperature region, where dxs(Fl,T)/dT > 0. As 
the temperature is increased, the system gets into a re- 
gion for which dx B {H,T)/dT < with a FM fraction 
greater than the lower limit for xb, so it remains blocked 
without changes in the magnetic response, a fact char- 
acterized by the plateau observed in M(T) for the high 
temperature region. This non-equilibrium phase diagram 
was constructed for the particular measurement proce- 
dure employed for the acquisition of M(H) and M(T) 
data. A modification in this procedure (for instance, by 
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FIG. 5: (color online). State diagram of the system in the 
x — T plane, for H = 1.3T, resulting from Eq. (9). Re- 
gions labeled dynamic indicate that the system can evolve in 
the measured time. Those labeled frozen indicates that the 
system is blocked. Data obtained from M(H) measurements 
at different temperatures ( up triangles) and from M(T) at 
H=1.3 T after ZFC (down triangles) are also shown. The 
dashed lines show the new phase diagram boundaries if the 
measuring time is increased by one order of magnitude. 

changing the time spend at each measured point) will re- 
sult in a change in the factor v cxp , with the consequent 
change of the boundaries in the phase diagram. For ex- 
ample, the effect of increasing the measured time by one 
order of magnitude is shown by dashed lines in the phase 
diagram of Fig. 5. 



IV. CONCLUSIONS 

In conclusion, we presented a thermodynamic phe- 
nomenological model for a global description of the phase 
separated state of manganites. The construction starts 
with the calculation of the free energies of the homoge- 
neous FM and CO phases. The free energies obtained 
turned out to be very close in value: the difference was 
of the order of the magnetic energy for intermediate fields 
in the whole temperature range investigated. The phase 
separated state is introduced by considering a uniform 
spread of the free energy density of the FM phase, and 
the dynamic behavior is included within a scenario in 
which the evolution of the system is determined through 
a cooperative hierarchical dynamic with diverging energy 
barriers as the system approaches equilibrium. The main 
success of the model is to provide an understanding of the 
response of the phase separated state when both temper- 
ature and magnetic field are varied, being able to repro- 
duce the dynamic and static properties of the system un- 
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der study. The same methodology can be also ap plie d to 
other compounds sharing similar phase diagrams 
and properties, especially those displaying abrupt field 
induced transitions at low temperatures. \2dL l26l l27j The 
key factors to determine the free energies of the homo- 
geneous phases is the possibility to measure the specific 
heat of each phase separately, taking advantage of the 
existence of blocked states, and the measurement of the 
temperature reached by the compound under study after 
the CO-FM abrupt transition at low temperature. This 
last value and the field at which the abrupt transition 
occurs are the key parameters to determine the homoge- 
neous ground state of the system at zero applied field. 

This work was partially supported by Fundacion An- 
torchas (Argentina) and CAPES, CNPq, and FAPERJ 
(Brazil). 
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